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ABSTRACT 

The vicinity of the unidentified EGRET source 3EG J1420-6038 has under- 
gone extensive study in the search for counterparts, revealing the energetic young 
pulsar PSR J1420-6048 and its surrounding wind nebula as a likely candidate for 
at least part of the emission from this bright and extended gamma-ray source. 
We report on new Suzaku observations of PSR J1420-6048, along with analysis 
of archival XMM Newton data. The low background of Suzaku permits mapping 
of the extended X-ray nebula, indicating a tail stretching ~ 8' north of the pul- 
sar. The X-ray data, along with archival radio and VHE data, hint at a pulsar 
birthsite to the North, and yield insights into its evolution and the properties 
of the ambient medium. We further explore such properties by modeling the 
spectral energy distribution (SED) of the extended nebula. 

Subject headings: pulsars: individual (PSR J1420-6048) - X-rays: general 



1. Introduction 



The campaign to identify 3EG J1420-6038 has revealed sources across the electromag- 
netic spectrum from radio to VHE 7-rays. The co mplex of compact a nd extended radio 
sources in this region is referred to as the Kookaburra (iRoberts et al.lll999l ) , and covers nearly 
a square de gree along the Galactic plane. Within a Northeasterly excess in this complex la- 
beled "K3" iD'Amico et all fepOlh discovered PSR J1420-6048 (hereafter J1420), a young 
energetic pulsar with period 68 ms, characteristic age r c = 1 3 kyr, and spin down energy 
E = 1.0 x 10 37 ergs _1 . The NE2001 dispersion measure model (ICordes fc Lazioll2002l ) of this 



Roberts et al 



(I2001M ) 



pulsar places it 5.6 kpc distant. Subsequent ASCA obse rvations b y _ 
revealed extended X-ray emission around this pulsar, and Ng et al. ( 20051 ) further examined 
the K3 pulsar wind nebula (PWN) with Chandra and XMM-Newton, r esolving a bright 
i nner n ebula along with fainter emission extending ~ 2' from the pulsar. lAharonian et al. 
( ]2006bl ) report on the discovery of two bright VHE 7-ray sources coincident with the Kook- 



aburra complex. HESS J1420-607 is centered just north of J1420, with best fit extension 
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overlapping the pulsar position. The other H.E.S.S. source appears to corr espond to the 



Rabb it nebula h alf a degree southwes t, which is also observed in the radio ( [Roberts et al. 



19991 ) and X-ray flRoberts et al.ll2001ah. Most recently. PSR J1420-6048 was detected by the 



Fermi Large Area Telescope (LAT) ( lAbdo et al. II2009I ). This crowded region clearly merits 
further study, and we report on new X-ray results obtained with Suzaku and XMM-Newton, 
as well as SED modeling of the K3 nebula. 



2. Data Analysis 

The Suzaku pointing (obsID 503110010) occurred on January 11-12 2009 for a total of 
50.3 ks. We utilize the standard pipeline screened events, and analyze the XIS front side 
(XIS0 and XIS3) and back side (XIS1) illuminated chips with XSelect version 2.4. We also 
obtained recent archival XMM data to augment the Suzaku data; observation 0505840101 
occurred on February 15 2008, for 35.0 ks, while observation 0505840201 added 5.6 ks. The 
second data set has a slightly different CCD placement, and suffers from high background, 
so we only use the 35.0 ks of data. We apply the standard data processing, utilizing SAS 
version 9.0. After screening the data for periods of high background 19.9 ks remain with the 
MOS chips. The PN chip suffers greatly from flaring, and we discard this data. Spectral fits 
are accomplished with XSPEC version 12.5. 



2.1. Broadband Morphology and Point Sources 

Suzaku X-ray emission is peaked in the vicinity of the pulsar, with a bright halo ex- 
tending ~ 3' and a fainter tail extending north ~ 8'. A number of other excesses of emission 
correspond to point sources, as discussed below. Figure 1 shows the Suzaku data in the 
2-10 keV band, which highlights the extended PWN emission to the north. Also depicted is 
the XMM exposure, clearly showing a number of point sources, though no obvious extended 
emission is apparent. 

To identify X-ray point sources we use the SAS source detection function edetect_chain 
on the XMM MOS chips and search in two energy bands of 0.3-2 keV and 2-10 keV for 
sources with a probability P < 10~ 13 of arising from random Poissonian fluctuations. The 
source detection algorithm also attempts to determine source extension via a Gaussian model, 
though all detections are consistent with a point source. Counts are therefore extracted from 
a 15 pixel ( 1 6.5") radius circle. A total of 8 sources pass this test, 4 of which also appear 



m 



Ng et al.l ( 120051 ): PSR J142-6048 (source 5 in our dataset), the X-ray sources denoted 
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Fig. 1.— Left: Suzaku XIS0+XIS1+XIS3 2-10 keV image overlaid with H.E.S.S contours 
and smoothed with a 40" Gaussian. PSR J1420-6048 is indicated by the white cross. Right: 
XMM EPIC 0.2-12 keV mosaic showing named spectral extraction regions and numbered 
point sources, smoothed with a 10" Gaussian; no extended tail emission is apparent. 



star 1 (source 1) and star 3 (sou rce 2), and another point source to the southeast (source 



3) unlabeled by (jNg et al.l 120051 ) but visible in their XMM exposure. This source to the 
southeast is also a field star, as it appears quite bright in DSS2 red images. Of the four 
remaining sources, only one, a hard bright source 8.5' north of the J1420 labeled source 7, 
lacks an optical counterpart. Source 7 also overlaps a radio hotspot to the north. Below we 
list the properties of these sources, defining the hardness ratio as: HR=(Chi — Ci G )/ (Chi + Ci ) 
where Ci Q and Chi are MOS counts in the 0.3-2 keV and 2-10 keV bands, respectively. It is 
worth noting that PSR J1420-6048 is only the fifth brightest point source in the XMM field 
and that all 8 XMM sources appear as excesses in the soft band Suzaku data as well. All 
point sources are quite soft, save for J 1420 and source 7. 

On a larger scale, extended emission is observed in all wavebands. Australia Compact 
Telescope Array (ATCA) observations within the error ellipse of 3EG J142 0-6038 (which 



i s bro ad enough to encompass both the K3 wing and the Rabbit nebula) by iRoberts et al. 



(119991 ) revealed the "K3" excess, a resolved knot of emission surrounding the pulsar of flux 
density 20 mJy at 20 cm with index a = —0.4 ± 0.5. Adjacent is the "K2 wing," with 1 
Jy at 20 cm and index of —0.2 ± 0.2. Closer inspection of the both the 13 cm and 20 cm 
continuum maps reveal that J 1420 lies on the southeastern rim of an apparent radio shell ~ 3' 
in radius. This shell is also apparent in the SUMSS 843 MHz map of the region. The center 
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Table 1: XMM Source Properties 



No (* ) 




R.A. 




Dec. 




Pos. Err." 


Counts 


HR 


1 (Star 1) 


14 


19 : 


11.52 


-60 


: 49 : 


34.00 


0.26 


462 ± 27 


-0.99 ±0.026 


2 (Star 3) 


14 


19 : 


31.48 


-60 


: 46 : 


20.29 


0.39 


146 ± 18 


-0.85 ±0.12 


3 (Unlabeled) 


14 


20 : 


22.72 


-60 


: 53 : 


21.47 


0.45 


164 ± 19 


-1.00 ±0.069 


4 


14 


19 : 


17.61 


-60 


: 45 : 


23.45 


0.61 


123 ± 17 


-0.81 ±0.12 


5 (PSR J1420-6048) 


14 


20 : 


08.19 


-60 


: 48 : 


14.85 


0.63 


150 ± 20 


0.95 ±0.072 


6 


14 


20 : 


40.75 


-60 


: 41 : 


20.22 


0.79 


40 ± 11 


-0.79 ±0.30 


7 


14 


20 : 


09.78 


-60 


: 39 : 


42.86 


0.85 


62 ± 14 


0.79 ±0.18 


8 


14 


19 : 


35.85 


-60 


: 42 : 


11.39 


1.16 


34 ±11 


-0.68 ±0.36 



(|Ng et al.ll2005l ) counterpart 



of this shell coincides with a dearth of emission in Spitzer 8 /im maps. We place extended 
PWN upper limits in the radio by measuring the 843 MHz, 20 cm and 13 cm flux densities 
from the entire shell (which is significantly larger than the X-ray extension), finding 0.61, 
0.75 Janksy and 0.58 Jansky, respectively. We deem these flux densities to be upper limits 
since the poor spatial resolution of Suzaku prevents ascertaining how the radio and X-ray 
emitting regions relate. We also remeasure the K3 excess and find flux densities of 15 mJy, 
19 mJy and 17 mJy at 843 MHz, 20 cm and 13 cm, respectively, consistent with the result 
of iRoberts et all dl999h V At higher energies HESS J1420-607 shines at 13% of the Crab 
(lAharonian et al.ll2006bl ). with photon index of 2.2 and extent of 3.3' centered 2.6' north of 
the pulsar. The H.E.S.S. spectrum is extracted from a 9.6' circle to minimize contamination 
from HESS J1418-609 (The Rabbit) 33' to the southwest. 



2.2. Spectrum 

The XMM point source identification described above is valuable since the broad Suzaku 
PSF hinders the disentanglement of point sources from extended emission. We are accord- 
ingly free to define extended regions for spectral analysis while steering clear of the X-ray 
point sources; we label these "Inner," "Halo," and "Tail". All spectral extraction regions 
lie external to the Suzaku 1' half power radius of the 8 point sources (except the very faint 
and soft source 8, which slightly overlaps the "Tail"). The Northernmost emission is encom- 
passed by the "Tail" region which extends north of the pulsar from 3' ~ 8'. Suzaku does best 
with large regions, and accordingly we capture the pulsar and surrounding PWN emission 
with a circular region of radius 3' ("Halo"). To best isolate the pulsar and inner PWN we 
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Fig. 2. — Left: 20 cm ATCA image s howing the possible c avity northwest of PSR J1420- 
6048 (indicated by a white cross) after iRoberts et al.l (119991 ). Suzaku contours are overlaid, 
with the northern blob due to source 7. Right: SUMMS 843 MHz image with H.E.S.S. 
contours overlaid. The bright emission to the south is the HII region G313.5+0.2. 



also use a 1.7' circle we call "Inner," which is the minimum size recommended for Suzaku 
spectral analysis. The Suzaku background region occupies the northern and eastern edges 
of the field of view. 

We extract Suzaku XIS spectra with XSelect from all three active XIS chips, while re- 
sponse files (both ARF and RMF) are generated with the xisresp script. The ASCA FTOOL 
addascaspec is utilized to combine the spectra and responses of the XIS front side illuminated 
chips (XISO and XIS3). The XIS1 back side illuminated chip possesses a markedly different 
response function from the front side illuminated chips, and so is analyzed in parallel rather 
than added to the other two chips. Finally, spectra are binned to a minimum of 30 counts 
per bin. The high XMM background prevents an adequate fit to the extended emission (flux 
errors of ~ 100%). We are able, however, to extract the spectrum from the "Inner" region 
encompassing the pulsar. To minimize calibration errors due to differing chips and distance 
from the pointing axis, we select an XMM background region on the northeast corner of 
the central chip containing PSR J1420-6048. We extract XMM spectral data with the SAS 
function evselect, create responses with the functions rmfgen and arfgen, and group to 30 
counts per bin. 



To best probe spectral index variations we make a global estimate for the best fit 
hydrogen column density. To maximize the signal to noise, we simultaneously fit nu for the 
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Suzaku "Halo" spectrum and the "Inner" XMM spectrum. For an absorbed power law model 
we find a n# of 4.1!q 4 x 10 22 cm -2 (90 % single parameter error). An individual fit to the inner 
PWN region indicates that an absorbed power law model provides the best fit to the data; 
addition of a blackbody, neutron star atmosphere, or thermal plasma (mekal) component 
does not improve the fit. Below 2 keV the fit to all regions significantly underestimates 
the flux. This feature is not well fit with a thermal plasma, neutron star atmosphere, or 
blackbody, and might simply constitute excess soft emission from the myriad faint sources 



in this crowded reg ion. The 4.1 



+0 j x 10 



22 



cm 2 column is marginally consistent with the 



value measured by iRoberts et all fl200fri ) of 2.2 ± 0.7 x 10 22 cm~ 2 using ASCAGIS data 
(3' radius), and matches well with the short Chandra ACIS exposure fit by iNg et al.l (120051 ) 
(2' radius aperture) of 5.4l 2,2 x 10 22 cm -2 (la err or). The nominal tot al Galactic H column 
i n this direction is estima ted as 1.6 x 10 22 cm -2 fjKalberla et al.l 120051 ) and 2.1 x 10 22 cm -2 
rtPickev fe Lockmanllll)9~oT ). 

With uh fixed at 4.1 x 10 22 cm -2 , we extract spectra from the three regions described 
above. Table 1 lists the results of fitting an absorbed power law to these regions in the 2-10 
keV range; errors are 90% single parameter values. We quote 2-10 keV fluxes to compare 
to previous work, minimize sensitivity to the hydrogen column density fit, and mitigate soft 
X-ray contamination in this crowded region. In the inner 100" region we fit simultaneously 
the XMM and Suzaku data, and find a power law index of 1.8 ± 0.1. The entire 3' "Halo" is 
fit with an index of 2.0 ± 0.1, implying a steeper power law near the edge of the nebula. To 
further explore this, we extract the spectrum from an annulus of radius 1.7' — 3' from J1420. 
For such annular regions the Suzaku response normalization proves unreliable, though the 
spectral shape is still of value. We find a power law index of 1.9 ± 0.1; this is indeed softer 
(though within errors) of the inner PWN. Further spectral softening is hinted at in the outer 
nebula "Tail " which boasts an index of 2.1 ± 0.1. 



Table 2: Spectral Fits 

Region Nn T abs. flux unabs. flux x 2 /dof 

xl0 22 cm- 2 / 2 _ 10 t / 2 _ 10 t 

Inner 4.1* 1.82 ±0.13 1.19±g;lo l-62±g;^ 333/345 

Halo 4.1* 2.00 ±0.10 2.45^jf 3.40+^ 219/309 

Tail 4.1* 2.14 ±0.16 1.26±gjg 1.7918;^ 138/224 



Simultaneous fit to Suzaku and XMM data 
*held fixed 

^2-lOkeV fluxes in units of 10~ 12 ergcm~ 2 s _1 



-7- 



Thes e spectral values are largely consistent (though with smaller errors) with previous 
data fits. iNg et al.l ( 120051 ) extract the spectrum of the inner nebula from a 2' circle, which 
gives (for riu fixed at 5.4 x 10 22 cm -2 ) Y = 2.3^gg (projected multidimensional 1 sigma 
error), unabsorbed 2-10 keV flux 1.3 ± 0.14 x 10~ 12 erg cm -2 s -1 (1 sigma single parameter 



error 
4.7 x 10 



Roberts et al.l (j2001bl ) measure a power law index of 1.6 ± 0.4 and a 2-10 keV flux of 



12 erg cm 2 s 1 for an extraction region of size 3'. 



3. SED Fitting 

Matching SED data points can help constrain physical properties of the source. To this 
end we we apply a two-zone time-dependent numerical model with constant injection lumi- 
nosity. We inject a power law spectrum of relativistic particles (either electrons or protons) 
with a high energy exponential cutoff into zone 1, and then evolve this spectrum over time 
according to radiation losses from synchrotron and IC (Klein-Nishina effects included). The 
resultant spectrum is subsequently injected into zone 2 and further evolved. Adiabatic cool- 
ing is ignored, given that offsetting heating effects of the SNR reverse shock interaction with 
the PWN may significantly alter par ticle energetics. The typical timescale for the reverse 



shock to reach the PWN is ~ 7 kyr ([Reynolds &: Chevalier! 1 1984 ). while the characteristic 



age of the pulsar is only 13 kyr, so this perturbation may still be felt. The complex nature 
of this interaction is also why we use a constant injection luminosity rather than allowing 
injection luminosity to vary with pulsar spin down power. Injection (and evolution) occurs 
in time steps much smaller than the assumed age. We consider a low density environment 
with density n = 1 cm -3 and photon fields comprised of CMBR, near IR, and far IR (though 
starlight is unimportant due to Klein-Nishina suppression of IC at high photon energies). 
Our two zones are defined as the inner 100" region surrounding the pulsar (zone 1) and the 
broader extended nebula (zone 2). Figures 3 and 4 show the SED of J1420. While we must 
select spectral extraction regions such that interference from stellar sources is minimized and 
region surface brightness remains significantly above background, the nebular flux extends 
beyond these regions. An estimate of the total nebular flux is essential if we wish to compare 
the X-ray flux with the VHE flux, which is extracted from a large ~ 10' circle. To estimate 
the total extended nebula flux, we extract 2-10 keV counts from a large ellipse 5' x 6' in 
radius which seems to encompass the majority of the X-ray flux (excluding point sources 
within this field) and find 40% more counts in this region than the combined "Halo" plus 
"Tail" minus "Inner" region. We therefore extrapolate that the total extended nebular flux 
is 40% greater than the combined "Halo" plus "Tail" minus "Inner" fluxes. Accordingly, for 
SED modeling of the extended nebula we plot the spectral points from the "Tail" region, 
scaled up to represent this estimated 2-10 keV flux of 5.0 x 10~ 12 ergcm -2 s _1 . 
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wn Figures 3 and 4 We plot the gamma-ray flux from HESS J 1420-607, along with radio 
and X-ray data. Also indicated are the Fermi LAT one-year and ten-year 5cr flux limit for a 
point source residing in a background 10 x greater than the extragalactic background, which 
is reasonable for the Galactic plane. 



3.0.1. Leptonic gamma-ray emission model 

For an X-ray spectral index approximately flat on the SED plot one cannot easily match 
the X-ray data to the H.E.S.S. data. A single electron injection component works, but only 
after significant fine tuning of injection and ambient medium parameters. A dual component 
injection comprises the other alternative, though the attentive reader might note that the 
northern X-ray nebula spatially corresponds quite well to the H.E.S.S. position and size, 
hinting at a single electron population responsible for both types of emission. If one assumes 
all particles evolve in the same spatial region (and hence same magnetic field and ambient 
photon field) one cannot independently fit the X-ray data to one component and VHE 
data to another significantly different component. Therefore, we refrain from displaying the 
unsatisfactory dual component fits, and focus on a single component model. 

A single electron component demands that the H.E.S.S. detection arise from IC scatter- 
ing of far IR photons, as the CMB provides insufficient seed photons to account for the VHE 
peak given the constraints imposed by the X-ray data. For the Galactic radius of J 1420, fa r 



IR photons typically peak at T = 25 K with a density ~ 1 eV cm -3 (IPorter et al.l 120061 ) . 
We find a better fit with a denser far IR photon field (2 eVcm -3 ), which is reasonable given 
the large scatter in ambient photon fields throughout the Galaxy. The starlight photon field 
is taken to have a temperature of 3500 K with density 1 eVcm -3 . 

The innermost 100" region is populated by young electrons emitting in a relatively 
strong magnetic field. In the ideal case of a magnetic dipole, the pulsar surface magnetic 
field can be estimated as B = 3.2 x 10 19 (P P) 1/2 G where P is in seconds. This gives a 
surface magnetic field of 2.4 x 10 12 G for PSR J1420-6048. If one assumes the magnetic 
field is dipolar out to the light cylinder, and then falls off as the inverse of radius past the 
light cylinder, the volume averaged magnetic field out to 50" from the pulsar is calculated as 
11 fiG. Without knowledge of where the bulk of the synchrotron emission stems from (given 
the poor spatial resolution of Suzaku) we therefore adopt a mean magnetic field of of 12 //G 
for the innermost region. Electrons responsible for the hard X-rays can diffuse out of this 
region in as little as 400 years, as shall be seen in Section 4. We therefore evolve the inner 
zone 1 electrons over 500 years. A power-law distribution of electrons with the standard 
index of 2 is injected, with a high energy cutoff at 400 TeV. For these parameters, we require 
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an energy in the innermost region of 9 x 10 45 d\ 6 erg in the form of leptons, where distance 
d = g?5.6 x 5.6 kpc. The electron index of 2 underpredicts the K3 radio flux excess by a factor 
of ~ 4. A slightly softer electron index allows us to fit the radio points, but increases the 
energy requirements by a factor of ~ 3 and results in an outer nebula radio flux butting 
up against the upper limits established above. Better constraints on the non-thermal radio 
emission from this region would aid greatly in differentiating between models. 

In a magnetic field of 12 /xG evolving over 500 years results in a lepton spectral break 
at ~ 100 TeV. This broken spectrum is subsequently injected into zone 2, the outer nebula. 
With this electron spectrum and for the photon field described above the outer nebula 
requires a total lepton energy of 9 x 10 i7 d§ s erg in order to fit the H.E.S.S. data. This 
energy requirement is reasonable for an extended P WN, and can be compared to the Vela X 



nebula possessing ~ 10 48 erg in the form of leptons (jde Jagerl 120071 ). Evolving the electrons 
over 12.5 kyr places the X-ray and VHE emitting electrons in the cooled regime for the 8 /xG 
field we select to match the X-ray data. 



3.0.2. Hybrid hadronic and leptonic gamma-ray emission model 
One can also model the VHE gamma-rays as the product of pion decay from proton 



proton interactions. The time scale for pion prod uction via p-p interactions is given by r x 



1.5 x 10 8 (n/lcm 3 ) 1 years (jBlumenthall Il970l ) ; this timescale is significantly greater than 



pp 



the expected age of the system, so the proton spectrum is treated as static. For an injected 
power law proton spectrum, we calculate t he photons from pr oton-proton interactions and 



subsequent 7r° and 77-meson decay following iKelner et al.l ( 120061 ). Proton-proton interactions 
also yield ^ mesons which decay into secondary electrons, which we evolve over 13 kyr. IC 
and synchrotron emission from the resultant secondary electron spectrum are subsequently 
calculated. 

Secondary electrons evolved over 13 kyr cannot possibly account for the radio or X-ray 
data from the extended nebula for plausible magnetic fields. We nevertheless indicate on 
Figure 4 the synchrotron radiation from secondary electrons. The X-ray data must therefore 
be attributed to synchrotron radiation from a primary electron population. We adopt the 
same two-zone lepton model described above, though with the outer nebula parameters 
tweaked slightly. We adopt the typical far IR field density of 1 eVcm -3 , reduce lepton 
injection energy by ~ 15 % in the outer nebula, and adopt a magnetic field of 9 /iG for 
this region. A magnetic field any greater than this is precluded both by the radio data 
and by increased synchrotron cooling, which steepens the synchrotron slope in the Suzaku 
range. With these parameters gammas from pion decay can account approximately equally 
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Log Energy (eV) 



Fig. 3. — SED for single electron component showing the radio, X-ray and VHE data. Solid 
lines indicate the fit to the broad extended nebula, while the dot-dashed line correspond 
to the inner PWN region. Thinner lines denote the individual IC components: starlight 
(dotted), CMB (long-dashed), and far IR (short-dashed). Flux in the Fermi LAT realm is 
significantly below the 1 and 10-year detection lines. 

with IC emission for the VHE data. The hadronic contribution to the H.E.S.S. detection 
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can be modeled by a a proton power law of index 2, cutoff at 200 TeV, and energy E = 
7 x 10 50 (n/1 cm -3 ) d 2 6 erg. The SED from this hybrid hadronic plus leptonic scenario is 
shown in Fig 4. We omit the inconsequential IC radiation from secondary electrons. 

4. Discussion 

X-ray flux extending 8' north of the pulsar indicates a spatial extent of R = 13 d 5 6 pc. 
Though this requires an average flow speed of under 1% the speed of light for a 13 kyr age, 
diffusion across magnetic field lines can be an extremely slow process. For a cooling limited 
electron source with Bohm diffusion, the diffusion timescale varies with size (6), magnetic 
field (B), and electron energy (E e ) as: 

tdiff ~ 68 (0/8') 2 dl 6 (B/lOfiG) (-E^/lOO TeV) -1 kyr (1) 

Therefore over the ~ 10 kyr lifetime of J1420, the inferred presence of distant VHE electrons 
is likely due to a different process than simple diffusion, as discussed in the next section. 

The observation of hard synchrotron X-rays requires an efficient particle accelerator, 
since in a transverse magnetic field the mean synchrotron photon energy (E y ) scales as: 

£ 7 « 2.2(£ e /100 TeV) 2 (£/10 AtG) keV (2) 

Such energetic electrons cool rapidly, however, on a timescale of: 

r sync « 820 (EJ 100 TeVy^B/ 10 fiCy 2 year (3) 

Therefore the cooling timescale for electrons radiating photons E 1 at mean energy is: 

T sync « 1200 (E 7 /lkeV)- 1/2 ( J B/10 / uG)- 3/2 year (4) 

Synchrotron cooling dominates over IC cooling for electron energies above 5 TeV for the 
external photon and magnetic fields chosen, so we consider only synchrotron cooling in the 
following discussion. The electron cooling break results in a similar break in the emitted 
photon spectrum, which scales with magnetic field and age (r) as: 

E^Break ~ 26 (t/13 kyr)~ 2 (I?/10 fiG)~ 3 eV (5) 

This break is apparent in the models plotted in Figures 3 (break at ~ 50 eV) and 4 (break 
at ~ 30 eV). 

The break location is also of key importance to IC modeling of the H.E.S.S. detection. 
Electrons upscatter blackbody photons (at temperature T) to a mean energy of: 

£ 7 « 2.2 (EJ10 TeV) 2 (T/25 K) TeV (6) 
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Fig. A.— SED for hadronic plus leptonic model. Dot-dashed lines indicate pion decay 
and synchrotron from secondary electrons for a density of n = 1 cm -3 , solid lines denote 
synchrotron and IC from the primary electron population. The long dashed line denotes 
synchrotron flux from the inner nebula, and the solid line intersecting the H.E.S.S. points 
marks the total extended nebula flux. In this scenario Fermi detection could be expected 
after a few years. 
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where we scale T to the value adopted for the far IR field. HESS J1420-607 appears to peak 
in flux at w 0.5 — 2.5 TeV; we expect this peak to result from a cooling break rather than an 
exponential cutoff in the electron spectrum due to the need for significant numbers of VHE 
electrons to synchrotron radiate in the X-ray regime. Inserting the value of E 1 = 0.5 — 2.5 
TeV into Equation 6 we see that at the break: 

E e , Break « (5 - 11) (T/25 K)- 1 ' 2 TeV (7) 

This VHE break is primarily due to synchrotron cooling, so substituting Equation 7 into 
Equation 3 yields an estimate for the mean magnetic field strength in the nebula: 

B«(8- 11) (T/25 K)" 1/4 (r/13 kyr)- 1/2 /iG (8) 

While this range of magnetic field (which ignores IC cooling) should not be taken too seri- 
ously, the fact that the 8 /iG field we selected for single component modeling also matches 
the X-ray data lends some credence to this estimate. 

Comparing the merits of each model, we find that the leptonic model provides a slightly 
better fit to the data, without the cost of the additional spectral components of the hybrid 
model. In addition, in order for pion decay to account for the VHE emission completely the 
energy in protons exceeds 10 51 erg for a typical density of n — 1 cm -3 , and the fit to the lower 
energy H.E.S.S. data points is unimpressive. Even with > 10 51 erg in hadrons, a primary 
electron population required to match the X-ray data contributes significantly to the VHE 
flux via IC, as seen in Figure 4. Only by raising the ambient magnetic field to > 20//G and 
lowering the total energy budget of this electron population can the IC from these electrons 
be rendered insignificant. Yet for a magnetic field of this magnitude synchrotron cooling 
results in a 2-10 keV photon spectrum far steeper than our measured data, plus such a high 
magnetic field over such a large region would be surprising. The high energy requirement 
in hadrons, as well as the need for two, rather than one, highly tuned injection components 
makes the hybrid model a less appealing explanation for the observed extended nebula fluxes 
than a single electron population. In addition, further support for the leptonic model stems 
from its ability to account for the innermost X-ray flux with a younger population of similar 
electrons radiating in the stronger magnetic field near the pulsar. 

Early detection by the Fermi LAT of the K3 nebula would undermine this assumption 
of a single electron component, however, as in this model very little flux is expected in the 
GeV range. Furthermore, a photon index of a = (p + l)/2 (where p is the electron index) 
is expected, which for our injection would imply a hard LAT index of 1.5. Detection by the 
LAT is much more likely in the hadronic scenario, and one would expect a GeV index of 
k 2. A LAT index steeper than 2 cannot be matched by either mo del, and would imply 



the existence of another particle population, likely electrons. Indeed, Ide Jager et al.l (120081 ) 
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suggested the Vela X nebula could be modeled with two populations of electrons. Such a 
situation for the K3 nebula would allow for detection by the LAT, and might also explain 
the radio excess. LAT detection of the K3 nebula would also greatly aid in distinguishing 
between scenarios, though even LAT flux upper limits might preclude the hybrid model. 



4.1. Birthsite 

The radio, X-ray, and VHE data all hint at a pulsar birthsite to the north. The possible 
radio filament extending to the hotspot 8.5' to the North might mark the synchrotron trail 
of J1420's proper motion. The nature of the X-ray (source 7) and radio point source at 
the end of this tail remains mysterious, though an AGN provides a likely candidate. The 
radio "filament" might also mark the eastern boundary of the apparent radio shell mentioned 
in Section 2. Under this assumption, the pulsar could have originated 3' northwest at the 
center of this shell. With this assumption, electrons spewed off the pulsar and left behind 
by its motion should also synchrotron radiate in the X-rays and emit VHE photons via IC 
scattering off field photons. The observed X-ray tail to the North, and the H.E.S.S. position 
offset ~ 3' in this direction nicely support this hypothesis. If the pulsar indeed originated 3' 
distant this implies a reasonable velocity of 

v = 370 (9/3') d 5 . G 0/13 kyr)" 2 km s" 1 (9) 

The northern tail could alternatively indicate an asymmetric reverse shock returning from 
the South has crushed and displaced the K3 PWN, much like the situation observed in Vela 



X (IBlondin et al.ll200ll ). The pulsa r birthsite would remain unco nstrained by our X-ray ob 



servations in this scenario. Indeed, iJohnston fc Weisbergl (120061 ) measured the polarization 



angle of the polarization sweep, which they argue correlates with the velocity. The perpen- 
dicular of this line (marking the pulsar spin and likely proper motion axis if it is emitting 
in the orthogonal mode) points southwest back to the shell body of the Kookaburra ~ 11' 
away. Association would require a very high space velocity of 1400 4.6 (t/13kyr) _1 kms -1 . 



5. Conclusions 

The Suzaku data provides the first deep mapping of the asymmetric extended X-ray 
emission surrounding J1420. The spectral index appears to soften in the outer nebula, as 
expected in a leptonic model with synchrotron cooling of electrons. SED studies of this 
nebula allow for a hybrid leptonic + hadronic model, yet favor a pure leptonic model if one 
takes into account energy requirements and the degree of tuning required to fit the hybrid 
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model. The outer nebula also yields information about the pulsar birthsite, hinting at an 
origin to to the North, though admittedly the northern X-ray tail could be due to reverse 
shock interaction with the PWN and therefore not correspond to the pulsar velocity. Further 
investigation of the apparent radio shell along with deeper Chandra imaging of the compact 
nebula might also help elucidate the origin of PSR J1420-6048. 

Acknowledgments: We thank Stefan Funk for useful comments on this paper. 
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